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ABSTRACT 

Herbig-Haro jets often show some degree of curvature along their path, in 
many cases produced by the ram pressure of a side-wind. We present simula- 
tions of both laboratory and astrophysical curved jets and experimental results 
from laboratory experiments. We discuss the properties and similarities of the 
laboratory and astrophysical flow, which show the formation of internal shocks 
and working surfaces. In particular the results illustrate how the break-up of the 
bow-shock and clumps in the flow are produced without invoking jet variabil- 
ity; we also discuss how jet rotation reduces the growth of the Rayleigh- Taylor 
instability in curved jets. 



Subject headings: ISM: jets and outflows — ISM: Herbig-Haro objects 
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Introduction 



Parsec-scale Herbig-Haro (HH) jets are one of the most spectacular phenomena con- 
nected to star formation. Observations of line emission features reveal an intricate succes- 
sion of bow shocks and knots along the jet, and a comp lex and rich dynamics ari sing from 
the interaction with the circumstellar environment (see iReipurth fc Bally I (120011 ) for a re- 
view). In the reflection nebula NGC 1333, in the Perseus molecular cloud, a number of 
bipo lar HH jets exhibit a di stinguishing C-shape morphology indicative of a steady bend- 
ing (IBally fc Reipurth Il200ll ). Less regular curvature is also observed in a number of other 
HH jets; for example in HH 30 a small side drift clos e to the jet source is followed fur- 
ther away by a sudden bending (lAnglada et al. 1 120071 ). In general, the curvature in jets 
has been linked either with the motion of the jet sources relative to the ambient medium 
or with the presence of a widespread outflow; both cases giving rise to an effective trans- 
verse ram pressure (cross-wind) which curves the jet. In fact, numerical simulations of jet 
bendi ng produced by ma gnetic fields seem to require field strength that are unrealistically 



high (IHurka et al. Ill999l ). For the jets in the Orion Nebula, it was also suggested that 



the bending may be caused by th e 'rocket' effect induced by the photo-io nizing UV radia- 



tion fr om nearby high-mass stars (IBally fc Reipurth 1120011 ); recent work by lKajdic fc Raga 



(120071 ) has investigated the influence of an ionizing photon flux on the jet-wind collision in 
the HH555 jet. The presence of curved jets is also common to other environments, such as 
in active galactic nuclei where they are observed as narrow- angle-tail radio sources (O'Dea 
and Owen 1986), or the pulsar jets in the Vela nebula, which are curved under the com- 
bined action of the wind within the supernovae remnant and the proper motion of the pulsar 
(Pavlov, Teter et al. 2 003). Under simplifying assumption the curv ature of the jets can be 
described analytically (IBegelman et al. Ill979l ; ICanto &: Raga Ill995l ). however these models 
cannot detail the shock structures and the intricate dynamics developing in the interaction, 
and that were observed in a number of numerical studies. For the narrow- angle-tail sources 
(IBalsara fc Norman Ill992l ) internal shocks were shown to be responsible for the formation 
of density enhancements (knots) inside the jet beam which was then disrupted through the 
combined effect of Kelvin-H elmholtz and Raylei gh- Taylor instabilities. HH jets were first 
investi gated numerically b y iLim fc Raga I (119981 ) for different incident angles of the wind 
and by iMasciadri fc Raga I (120011 ) for the case of irr adiated jets. More recently the jet-wind 



collis ions have also been studied in the laboratory (ILebedev et al. 1120041 ; lAmpleford et al 
20071 ) with supersonic (Mach number ~ 20 — 40) jets interacting with a fast, radiatively 



ablated plasma that acted as a cross-wind. The laboratory results clearly demonstrated a 
steady curvature in the jets under the impinging cross-wind and the appearance of internal 
shocks. 



In the present work we combine consistently the laboratory modelling (experiments and 
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simulations) with simulations of curved HH jets, drawing a parallel between the observed 
dynamics in the two systems and, with in mind future experiments, we also discuss the 
effects of rotation on curved jets. 



2. Laboratory Jets 



Laboratory jets are produced on the MAGPIE pulsed-power generator ([Mitchell et al 



19961 ) using a conical cage of micron-sized metallic wires (Al, Fe, W, etc) whi c h is s ubject 
to a 1 Mega Ampere current, with a rise-time of 250 ns (see iLebedev et al. I (120051 ) for a 
review). Figure 1 shows a schematic of the experimental set-up; the ohmically heated wires 
rapidly vaporize and turn into streams of plasma accelerated by the jxB^ force towards the 
axis of the array. The collision of the streams forms a standing conical shock, which redirects 
the plasma axially into a jet. Radiative cooling rapidly reduces the temperature in the jet 
which attains high Mach numbers (M~ 20 — 40) and small opening angles. The interaction 
of the jet with a side-wind is replicated in the laboratory by placing a plastic foil in the jet 
propagation region. XUV radiation from the standing shock radiatively ablates the foil and 
generates a fast moving plasma (~ 30 — 50 km s _1 ), providing the ram pressure to bend 
the jet. By changing the position, angle and length of the foil the characteristic parame- 
ters of the impinging wind can be altered, an d different regimes of the interaction studied 



(ILebedev et al. I l2004j ; lAmpleford et al. 1 120071 ). Similar to the HH jets, the characteristic 



cooling time in the laboratory jets is much smaller than the hydrodynamic time, this ratio 
is > 0.1, and radiation losses play a significant role in the energy balance of the jets. The 
numer ical simulations are performed with the thre e-dimensional resistive MHD code GOR- 
GON (IChittenden et al. I I2004J ; ICiardi et al. 1120071 ). the ion and electron energies are solved 
separately and include electron and ion thermal conductions. The coupling of the two energy 
equations due to collisional heat exchange is included as a sink/source term. The plasma 
is assumed to be in local thermodynamic equilibrium (LTE) and the average ionization is 
calculated by a Thomas-Fermi model. The plasma is optically thin and a radiation losses 
term is included in the electron energy equation. 

Two-dimensional axis-symmetric MHD simulations of the whole conical wire array were 
performed to obtain the mass and momentum fluxes of the jet. The calculated time and 
spatially varying profiles are then used as boundary conditions to inject the jet into the three- 
dimensional (3D) computational domain. Typical axial jet velocities and ion densities vary 



in th e range Vj ~ 100 — 200 km s and rij ~ 10 — 10 cm respectively (ICiardi et al 



20021 ) ; the jet is injected with an initial temperature Tj = 15 eV ~ 175 x 10 3 K. The 



material of the wires, and thus the jet, is tungsten. In the jet and interaction region the 
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magnetic field is negligible and in the 3D simulations presented here only the gasdynamic 
equations are integrated in time. The physical conditions of the cross-wind, produced by 
photo-ablation of the foil, depend on the incident radiation flux coming from the conical 
shock and are also time-dependent. However the simulations show th at within the range of 



the expected and measured wind parameters (lAmpleford et al. 1 120071 ). the dynamics of the 



interaction remains substantially the same and a detailed treatment of the time-variation 
of the wind is not essential. For the present simulations the injected wind has a constant 
velocity (v w = 60 km s _1 ) and temperature (Tj = 30 eV), the characteristic ratio of the jet 
velocity (which we remind is time-dependent) to the wind velocity Vj/v w ~ 2 — 4. Two wind 
density profiles are discussed, one is constant n w = 1.2 x 10 19 cm -3 (case A) and the other 
is exponentially increasing (case B). However for both cases the density contrast ratio is in 
the range rij/n w ~ 0.1 — 10; the wind is supersonic and pressure gradients play no significant 
role in bending the jet. 

Simulation results for the laboratory jet- wind interaction are shown in Figure 2 as XUV 
emission maps (photon energies hv > 30 eV) integrated along the line of sight. The emit- 
ted radiation includes bound-bound, free-bound and free-free transition and it is calculated 



assum ing LTE within a screened hydrogenic model (IStehle fc Chieze 1 120021 ; iMichaut et al. 



20041 ) . Self-absorption of radiation and scattering are not taken into consideration. In the 
images the jet is injected through the lower boundary and the wind is injected on the left 
boundary. The collision of the wind with the jet generates a bow shock which envelopes the 
whole jet. As the bow shock is advected downstream with the wind, it develops a highly 
asymmetric shape, with the upwind side of the jet cocoon effectively disappearing. An 
oblique shock forms in the jet body which begins to bend and redirect the jet flow sideways. 
It is clear that the jet bends considerably more in case A where the time-averaged ram 
pressure of the wind is about an order magnitude higher than in case B. As the jet curves, 
momentum transferred from the jet to the working surface rapidly decreases, until the work- 



ing surface effectively detaches from the jet; a result described by (IMasciadri fc Raga 1120011 ) 
in the context of simulations of the jet HH 505. The subsequent propagation of the working 
surface is ballistic and as it entrains more ambient mass it slows down and its emission rapidly 
decays. The presence of internal shocks in the jet beam is visible in both the simulations and 
experiments. The experimental XUV self-emission image (Figure 3) shows the formation of 
a new working surface half-way up the image and enhanced emission inside the jet beam. 
These shocks are produced by perturbations developing along the jet-wind interface and 
their effect is to slow down and redirect the jet sideway, allowing the un-shocked flow ahead 
to detach. The effect is particularly pronounced for the case, considered here, where the jet 
momentum increases in time. This produces a narrowing of the angle of the main oblique 
shock (with respect to the z-axis) and causes the jet to bore a new, straighter channel in the 
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ambient medium. Also a time- varying wind would have the same effect. Nevertheless the 
formation of new working surfaces was seen in simulations of uniform laboratory jets and 
wind, and in the simulations of astrophysical jets presented next. 



3. Jets from young stars 



The relevance of laboratory experiments to astrophysical jets, and for the case dis- 
cussed here the jet-wind interaction, rests on the ability of producing both an adequate 
representation of the dynamics of the astrophysical system and, m ore importantly obta ining 



dimensionless parameters in the appropriate astrophys ical range (IRvutov et al. Ill999l ). For 



the present experiments such conditions are well met (ILebedev et al. I l2005h and we expect 
much of the dynamics observed in the laboratory to be valid for astrophysical jets. 

Expected wind velocities vary from a few km s _1 for the jet-wind interaction associated 
with rel ative motions of TTauri stars with respect to the surrounding environment, (see for 
example Jones fc Herbig ( 19791 )) to typically higher velocities fo r irrad iated jets: best fits 
to HH505 Ha emission maps we re obtained by iMasciadri fc Raga I (120011 ) for a wind velocity 
of 15 km s _1 and estimates by iBally &: Reipurth I (120011 ) give wind velocities in the Orion 
nebula and in NGC1333 of ~ 10 — 20 km s _1 . Here we present simulations of HH jets in a 
cross-wind but do not model a specific jet. Instead we maintain the ratios of jet and wind 
velocities, and densities in a similar range as those obtained in the laboratory system. Specif- 
ically we take for the jet Vj = 100 km s -1 , rij = 1000 cm -3 and for the wind v w = 25 km 
s -1 , n w = 100 cm -3 . Although these values are consistent with observations, the somewhat 
larger wind velocity was chosen to produce enough bending in the jet while at the same 
time maintaining as much as possible of the interesting dynamics within the computational 
domain. Lower v w produce qualitatively similar results but because of the smaller curvature 
radius require considerably more computational resources. The initial jet radius rj = 50 
AU and the resolution is 4 AU over the entire Cartesian grid. The initial temperature in 
both jet and wind is 5000 K. The HH jet simulations are performed with an appropriately 
modified version of our laboratory code, where we now follow the time-dependent ioniza- 
tion of hydrogen. We take into account the recomb i nation and collisional ionization and 
use the rate coefficients as tabulated in iRaga et al. I (120071 ) (and references therein). For 
temperatures above 15000 K cooling is impleme nted by a function appropriate for interstel- 



lar gas composition (IDalgarno fc McCray 



19721 ). However for temperatures below 15000 K 



cooling is calculated by including collisional excitation, collisional ionization and radiative 
recombination of hydrogen, and the collisional excitation of O I and O II. Because of charge 
exchange the neutral and singly ionized populations of oxygen are assumed to follow closely 
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those of hydrogen (IHartigan &: Raymond Ill993l ); the atomic abundance of oxygen is taken 



to be 5 x 10 4 that of hydrogen 

The HH jet evolution is shown in Figure 4 in a time series of column density maps 
spanning 500 years of evolution. In general the overall jet dynamics is similar to that seen in 
laboratory jets. The presence of internal oblique shocks in the jet is apparent both in Figure 
4 and also in Figure 5c which plots the Ha emission. The simulations illustrate the formation 
of a number of bow-shaped shocks and knots in the resulting flow. It is worth noting that 
the dumpiness seen in the jet and the break-up of the bow shock arise entirely from the 
interaction with the cross-wind and are not imposed by means of a variable jet injection 
velocity. The latter is often used in numerical modeling to reproduce the knotty emissions 
features observed in the jets and which are thought to be caused by small velocity variations 
in the flow. The present simulations and more importantly the experiments, suggest that 
these internal shocks can also be produced by a wind impacting on a jet and are a dynamical 
observable feature dependent on the environment and not on the jet velocity profile. In 
addition to a steady cross-wind, the relative motion of the jet with respect to the interstellar 
medium, such as due to jet precession or the localized presence of transverse winds, may be 
also responsible for the formation of some of the knots observed. The detachment of the 
main working surface, as the jet body is increasingly curved under the action of the wind, 
occurs during the first 200 years and it is particularly evident in Figure 5b, which shows a 
passive tracer of the jet material at 370 years. The contact discontinuity between the jet and 
wind plasmas clearly shows perturbations along its surface which are ultimately responsible 
for producing the internal shocks in the jet. It is also evident that the jet is 'naked': the 
upstream side of the jet is ess entially in contact with th e wind and the cocoon is absent. 



As pointed out in the work of [Balsara fc Norman I (119921 ) the jet is liable to the growth of 



the Rayleigh- Taylor instability (RTI). We find that at the interface between the jet and the 
shocked wind, density and pressure gradients are in opposite directions, also the effective 
acceleration experienced by the jet material is in the direction of the centrifugal force and it 
is pointing into the density gradient, making the interface unstable. The development of the 
RTI in the jet beam is shown in Figures 5e and 5d which plot slices across the computational 
box (x~y plane) of the atomic number density. The wind is coming from the left and the jet 
is traveling out of the page. The arrows on the left-most panel indicate the height where 
the slices were taken. A well defined bow shock around the jet is clearly visible and so is 
the shock internal to the jet on the upwind side (Figure 5e). This is where the RTI can be 
seen to grow. The initial perturbation develops the characteristic bubble and spike structure 
(Figure 5d), with the tip of spikes rolling- up on the sides due to Kelvin- Helmholtz instability 
(KHI). The RTI tends to split the jet into well defined filaments; in addition as it penetrates 
deep into the jet the peak emission inside the jet shifts to the down-stream side, resulting 
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into a broadening of the emission across the whole diameter. However we might expect that 
any departures from uniform jet, wind and symmetry would modify the dynamics of the 



An example of particular interest is given by the simulation of a rotating jet shown in 
Figures 5f-j. The jet is taken to rotate as a solid body and we set the azimuthal velocity at 
its boundary = 0.5c s , where c s is the initial sound speed in the jet. The other parameters 
of the wind and jet are as those of the non-rotating case. The evolution is indeed similar, 
however the cross section of the rotating jet is somewhat larger due to the increased radial 
expansion and the density is correspondingly lower. The result is that the rotating jet tends 
to have a higher curvature (see for example Figures 5b and 5g) and a shorter propagation 
length. The most interesting difference though, can be seen in the slices across the com- 
putational domain (Figure 5i and 5j). Because of rotation (anti-clockwise in the images) 
perturbations driven by the RTI are advected around the surface of the jet and are rapidly 
sheared seeding the KHI; however for the instabilities are confined to a narrower region of the 
jet. The Ha emission maps in Figure 5c and 5h show that the internal shock for the rotating 
jet case remains narrower and the flow laminar for much longer. Thus rotation, combined 
with the action of the cross wind, reduces the growth of the RTI and the disruption of the 
jet. For the simulations presented here, the rotating jet survives up to z ~ 2000 AU as 
opposed to z ~ 1400 AU for the non-rotating case, above those heights the inabilities have 
grown to cover the whole cross sectional area of the jet. 



This work illustrates how combining laboratory experiments, numerical simulations of 
these experiments and the simulations of astrophysical systems can offer significant insights 
into the phenomena studied. In the case of the interaction of HH jets with a cross-wind, 
the results show the presence of oblique shocks and knots in curved jets, the formation of 
working surfaces, the detachment and fragmentation of the bow shock. In particular some 
of the observed features in HH jets may arise as a consequence of the curvature and the 
subsequent instabilities developing in the jet. However, it remains an open question what 
would be the minimum curvature required to produce them, especially for large values of the 
ratio PjVj/p w v w , where computations become very demanding. The simulations show the 
development of the RTI in curved jets and how it may be partially quenched by rotation, 
which rapidly shears the RT modes. To determine the observational signatures of such effect 
and for direct compar ison to observation, m ore realistic angular velocity profiles need to be 



RTI. 



4. 



Discussion 



used (see for example ILery &: Frank 
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A simple estimate of the RTI growth rate is obtained by assuming the plasma to be 
incompressible, its classical value is then given by to = \J Akg where k is t he wave number 



of th e perturbation and A = (p~j — p w ) / (p~j + p w ) is the Atwood number (I Chandr asekhar 



196ll ) evaluated with the local densities across the interface. For the jets we take the effective 
acceleration g on the jet fluid to be equal to the centrifugal acceleration: g = KVj, where 
k = 1/R and R is the radius of curvature. Thus for the simulated HH jets we have: Vj ~ 10 5 
m s" 1 , R ~ 3000 AU, k = 2n/rj ~ 8 x 10~ 13 m -1 , A ~ 0.3; where we have taken the 
wavelength of the perturbation to be the jet radius. The characteristic growth time is then 
t~rti = 7j ~ 70 years which compares very well with the growth time of the instability seen 
in the simulations. Laboratory jets also are susceptible to the RTI. However while it was 
seen in simulations, it is not clear if the filamentary structures observed in the experiments 
(cf. Figure 3) are caused by the growth of the RTI into the jet body. In general for the 
experimental jets Vj ~ 1.5 x 10 5 m s _1 , R ~ 0.05 m, k = 6300 m _1 , A ~ 1; giving a 
characteristic growth time, Trti ~ 150 ns, which is comparable to the propagation time of 
the jet in the cross-wind. 

A useful relation to estimate the growth rate of the RTI can be simply derived by using 
the analyti cal expression of the c urvature radius at the stagnation point for an isothermal 



jet given in ICanto fc Raga I (119951 ) . R = r^v 2 jV,J L c s 1 (pj/ p w ) 1 ^ 2 ■ The growth rate can then be 



3 w 

written as 

u = {2^A 1 l\ Pw /p J fl\c s v w f/\~ 1 

which also shows that within these simplifying approximations to does not depend on the jet 
velocity. The sound speed c s should be evaluated in the region of the jet body in pressure 
equilibrium with the shocked wind. Thus using reasonable values for the perturbed sound 
speed, c s ~ 5 — 15 km s -1 , we find for the simulated laboratory and astrophysical jets 
t rti ~ 150 — 250 ns and t rti ~ 100 — 200 years respectively; in good agreement with 
the previously estimated characteristic growth times. To observe the full development of the 
RTI in the laboratory jets we expect that experiments with longer interaction times (> trti) 
will be required. Such work is planned for the future and it will be coupled with the study 
of the effects of rotation on the jet. In fact, by using a mo dification of the conical a rray 



described in the present work we have recently demonstrated ( lAmpleford et al. Il2007bl ) the 
possibility of producing rotating jet in the laboratory with ratios of the axial to azimuthal 
velo city, vjvs ~ 5 — 10 , of the same order as those reported in recent observation of HH 



jets flCoffev et al. 1 120041 ) . 
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Fig. 1. — Schematic of the experimental set-up. 
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Fig. 2. — Synthetic XUV emission images of the laboratory jets shown at two times 290 ns 
and 320 ns from the start of the current pulse. Top two panels for case A and bottom two 
panels for case B (see text). The wind is injected through the left boundary. For case A 
the wind injection length is indicated; for case B the wind injection is over the entire left 
boundary. 
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Fig. 3. — Experimental time-resolved XUV image at 380 ns after the start of current. The 
foil in the experiments has an inclination of 30° with respect to the z-axis. 
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Fig. 4. — Time evolution of jet-wind interaction for astrophysical condition relevant to YSOs 
jets. The plots show on a linear scale the column density (g cm -2 ). The computational 
domain has a size 2004 x 4864 AU. 
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Fig. 5. — Comparison of the jet-wind collision for the non-rotating jet (panels a-e) and 
the rotating jet (panels f-j) cases at 370 years. Panels a-c and f-h show slices in the x-z 
plane of various quantities integrated along the y-direction; the size of the domain shown is 
2004 x 4864 AU. The combined jet-wind column density (g cm" 2 ) is shown in (a) and (f). 
Panels (b) and (g) show a tracer of the jet plasma on a logarithmic scale in arbitrary units. 
Panels (c) and (h) plot the Ha emission (erg cm" 2 s" 1 sr" 1 ) on a linear scale. Panels d, e, 
i and j are slices in the x-y plane across the computational domain and show the hydrogen 
number density (cm -3 ) on a logarithmic scale. The size of the domain shown is 2004 x 2004 
AU. The heights where the slices were taken is indicated by the arrows on the left-most 
panels; panels e and j correspond to z = 760 AU, panel d to z = 1261 AU and panel i to 
z = 1520 AU. 



